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The small scale dynamo and the amplification of magnetic fields in 
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ABSTRACT 

While our present standard model of cosmology yields no clear prediction for the initial mag- 
netic field strength, efficient dynamo action may compensate for initially weak seed fields via 
rapid amplification. In particular, the small-scale dynamo is expected to exponentially am- 
plify any weak magnetic field in the presence of turbulence. We explore whether this scenario 
is viable using cosmological magneto-hydrodynamics simulations modeling the formation of 
the first galaxies, which are expected to form in so-called atomic cooling halos with virial 
temperatures T v j r > 10 4 K. As previous calculations have shown that a high Jeans resolu- 
tion is needed to resolve turbulent structures and dynamo effects, our calculations employ 
resolutions of up to 128 cells per Jeans length. The presence of the dynamo can be clearly 
confirmed for resolutions of at least 64 cells per Jeans length, while saturation occurs at ap- 
proximate equipartition with turbulent energy. As a result of the large Reynolds numbers in 
primordial galaxies, we expect saturation to occur at early stages, implying magnetic field 
strengths of ~ 0.1 fiG at densities of 10 4 cm 4 . 
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1 INTRODUCTION 

Magnetic fields are ubiquitou s in the entire cosmos. They are ob- 
served in the local universe <Beck&Hoernej|l996l : lBeck||2007l : 
IChvzv & But a 2008: Keplev et al .1120101) and there is growing ev- 
idence for their existence at earlier cosmic times up to z ~ 4 
l lBernet et alJl2008l : lKronbere et al.ll20o3 : lMurphy||2009l) . Upcom- 
ing telescopes like LOFASQ] and SK^f will enable the detection 
of magnetic fields at higher redshifts. In present-day galaxies, the 
field strength ranges from a few to tens of pG. Their presence plays 
a vital role in the formation and evolution of many astrophysical 
systems. It is thus interesting to explore whether magnetic fields 
may have played a similar role at early cosmic times, particularly 
during the formation of the first galaxies. 

The first mag netic fields may have created via electroweak 
Bavm et aHll996t) or q uantum chromodynamics phase transitions 
Ouashnock et alj 19891) . Alternatively, they are generated by astro- 



physical mechanisms like the Biermann ba ttery effect or the Weibel 
instability during the structure form ation dBierman r] |l95d : IWeibeil 
ll959l : ISchlickeiser & Shuklall2003l) . Irrespective of their origin, the 
strength of magnetic fields predicted by these processes is many 
orders of magnitude lower than the present-day fields. Hence, they 
must have been amplified during cosmic history. 
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The simplest way to increase the magnetic field strength is 
gravitational compression under the constraint of flux-freezing. 
For a spherical collapse, it can be shown that B oc p 213 where 
B is the magnetic field strength and p is the gas mass den- 
sity. A more efficient mechanism is provided by astrophysical 
dynamos. In particular the small-scale dynam o can ampl i fy the 
small initial fields on very short timescales jKazantsevlll968l : 
iBrandenburg & Subram anian 20051). This process converts turbu- 
lent into magnetic energy, leading to magnetic field strengths 
in approximation equipartition with the turbulent energy. Nu- 
merous theoretical and numerical studies supp ort this idea in 
the context of cosmological structure formation ( Schleiche r et al.1 
20ld:lFederrath et al.l201 lllSur et alj20"l(j| : ISchober et al.l2012j|bl : 
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Turk et alj|2012h . vlodeling this process in numerical simulations 
is however challenging, as the turbulent eddies need to be well re- 
solved with at least 32 cells to provide a sufficiently large mag- 
netic Reynolds number Rm i Brandenburg & Subramanian! 120051 : 
ISuretalIl2oTol : iFederrath et alj 1201 ll) . The magnetic Reynolds 
number is given as Rm= VL/p, with V and L the characteristic 
turbulent velocity and length scale, and p. the magnetic diffusivity. 
The critical magnetic Reynolds number for dynamo action ranges 
from ~ 10 for Kolmogorov to more than 3000 for Burgers-type 
turbulence dSchoberetal]|2012iJ) . 

According to the current paradigm of structure formation, 
the first stars so-called Pop III stars, were form ed in mini-haloes 
with masses above > 10 s M Q at redshift 30 dAbel et alj 120021 : 
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iBromm & Larsonl2004l ; lYoshida et al.l2008l : IClark et al.l201 ll) . For 
these haloes, the efficiency of magnetic field amplification was 
demonstrated in various studies llSuret alJl20Tf] :l Federrath et al .1 
l201ll : lTurket alj|2012l ; lschober et alj|2012l) . The massive primor- 



dial haloes with virial temperature > 10 4 K were assembled around 
z = 15 -20 and hav e typical masses of > 10 7 M G dGreif et alj20O8t 
IBromm et aT1l2009l) . The main coolant in these massive haloes is 
the Lyman alpha radiation and they are known as atomic cool- 
ing haloes. The study of these haloes is of extreme interest as 
they are presumed to host the first protogalaxies which were the 
main drivers for the cosmic reionization as well as predecessors of 
present day galaxies. Exploring the magnetic field amplification in 
these haloes is more challenging, as a larger range of scales need 
to be bridged, and because of the higher Mach numbers involved 
in the problem. For such Mach numbers, a turbulent slope steeper 
than for Kolmogorov can be expected, implying a larger critical 
magnetic Reynolds number for the small-scale dynamo, and thus a 
higher minimal resolution. 

These atomic cooling haloes are considered as the poten- 
tial birthplaces for the formation of supermassive black holes 
at z > 6. The latter are presumed to power q u asars and have 
masse s of > 10 s M Q jFan et al.l 120031, 120061 : iMortlock et ail 
|20T 1). It is not yet fully comprehended how such objects could 
form in such a short span of time. Different mechanisms have 
been s uggest e d for the formation o f supermassive black holes 
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The essential condition for the direct collapse scenario is to 
avoid fragmentation. Molecular hydrogen is the only coolant in 
zero metallicity haloes which can bring the temperature of the 
gas down to a few hundred Kel vin and may induce fragmenta- 
tion by reducing the Jeans mass dLi et al.|[2003h . The existence of 
a strong UV Lyman-Werner flux photo-dissociates IT molecules, 
suppresses cooling and consequently may halt the fragmentation 
of a gas cloud. We consider a background UV flux with a ther- 
mal spectrum of T la( j = 10 5 K which quenches H2 formation via 
the Solomon process. These conditions create a suitable environ- 
ment for the growth of black holes. The presence of magnetic fields 
may favor this scenario by transferring angular momentum and ex- 
erting magnetic pressure depending on the ambient field strength. 
Hence, exploring the magnetization of atomic cooling haloes is in- 
timately associated with our understanding of the first galaxies and 
black hole formation. Indeed, strong rotation measures have been 
observed in quasars out to at least redshift 5, provi ding a further in- 
dicati on that magnetic fields are relevant early on dHammond et al .1 
l2012h . while detailed observations of the field structure in local 
AGN im ply significant con tributions to the transport of angular mo- 
mentum jBeck et al.ll2005l) . 

In this study, we aim to focus on the magnetization of massive 
primordial haloes, particularly the amplification of magnetic fields 
by the small scale dynamo. To accomplish this goal, we perform 
high resolution cosmological magnetohydrodynamical simulations 
with resolutions up to 128 cells per Jeans length. We also explore 



the saturation of magnetic fields due to the small scale dynamo 
for a Jeans resolution of 64 cells and examine its impact on the 
morphology of the halo. This study will enable us to comprehend 
the amplification of magnetic fields by the small scale dynamo in 
atomic cooling haloes and their potential role in the formation and 
growth of supermassive black holes. 

This article is organized in the following way. In the second 
section, we discuss the details of numerical techniques used in this 
study. In the third section, we discuss the results obtained. In sec- 
tion four, we explain the implications of magnetic fields in black 
hole formation. In the last section, we summarize the main findings 
and confer our conclusions. 



2 NUMERICAL TECHNIQUES AND SIMULATION 
SETUP 

The results presented in this work are obtained using the public 
version of the Enzo code. Enzo is an adaptive mesh refinement 
(AM R), parallel, grid-based cosmological magn etohydrodynamics 
code dO'Shea et al]l2004 iNorman et al.ll2007l) . The magnetohy- 
drodynamics (MHD) equations are solved using the Harten-Lax- 
van Leer(HLL) Riemann solver with piece-wise linear construc- 
tion. We enforce V.B = constraint employing the Dedner scheme 
dWang et alj|2008ll2010l) . 

The simulations are commenced with cosmological initial 
conditions generated from Gaussian random fields. We employ the 
inits package available with the public version of the Enzo code to 
create nested grid initial conditions. Our simulations start at red- 
shift z = 99 with a top grid resolution of 128 3 cells and we select 
the massive halo of 8 x 1 6 M Q at redshift 1 5 using the halo finder of 
iTurk et al. l|). Two initial nested levels of refinement are subse- 
quently added each with a resolution of 128 3 cells. Our simulation 
box has a cosmological size of 1 Mpc h -1 and is centered on the 
massive halo. In total, we initialize 6291456 particles to compute 
the evolution of the dark matter dynamics and have a final dark 
matter resolution of 300 M Q . The parameters for creating the initial 
conditions and the distribution of baryonic and dark matter com- 
ponen ts are taken from the WMAP seven years data djarosik et al .1 
l201ll) . We further allow additional 27 levels of refinement in the 
central 62 kpc region of the halo during the course of simulation. 
It gives us a total effective resolution of 0.17 AU at z = 11.5. The 
resolution criteria used in these simulations are based on the Jeans 
length, the gas over density and the particle mass resolution. The 
grid cells matching these requirement are marked for a refinement. 

The simulations conducted in this work mandated Jeans length 
resolutions of 16, 32, 64 and 128 (hereafter called J| 6 , J 32 , J64, J128) 
cells throughout their evolution. To include the primordial non- 
equilibrium chemistry, the rate equations of the following 9 species: 
H, H + , He, He + , He ++ , e", H", H 2 , H+ are self-consistently 
solved in the cosmological simulations. We make use of H 2 photo- 
dissociating background UV flux implemented in the Enzo code. 
An external UV field of constant strength 10 3 in units of J 2i is used 
in the simulations. We pre sume that such flux is generated from 
a nearby star forming halo ( Diikstra et alfeoolj) and is emitted by 
Pop III stars with a thermal spectrum of 10 5 K. We include several 
cooling and heating mechanisms like collisional ionization cool- 
ing, radiative recombination cooling, collisional excitation cool- 
ing, H 2 cooling as well as H 2 formation heati ng. The chemistr y 
solver used in this work is a modified version o f lAbel etal.ldl997l) ; 
lAnninos et al.l dl997l) . We stop our simulations after they reach the 
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Figure 1. This figure shows the radially binned spherically averaged radial profiles for the various Jeans resolutions as shown in the legend. The upper left 
panel shows the density radial profiles. The temperature radial profiles are depicted in the upper right panel. The bottom left panel shows the total energy radial 
profiles. The radial profiles for the total strength of magnetic fields are shown in the bottom right panel. 



maximum refinement level. The results at later stages would not be 
reliable. 



3 RESULTS 

In this section, we present the results obtained from cosmological 
MHD simulations fo r a halo of 3 x 10 7 M at redshift 1 1.5 selected 
from the sample of lLatif et al.l J2012h . The density perturbations 
begin to collapse under gravitational instability at a redshift around 
30. They merge with each other to form a halo in accordance with 
the hierarchical paradigm of structure formation. Gas falls into the 
dark matter potential and gets shock heated in the nonlinear phase 
of the collapse. Consequently, the temperature of the halo starts 
to increase. During the process of virialization, part of the gravita- 
tional energy is converted into thermal and kinetic random motions. 
This process continues until the halo reaches the state of virial equi- 
librium. 



3.1 Resolution study 

3.1.1 Radial profiles 

The radial profiles of the halo for different resolutions per Jeans 
length and the same central density are shown in figureQ] The max- 
imum density reached in the simulation is about 1CT 12 g cm 3 . The 



averaged density radial profile is shown in the top left panel of fig- 
ureQ] It can be seen that the density follows an R~ 2 behavior which 
is expected for close to isothermal collapse. The small bumps in 
the profiles for the J 32 and J 64 cases indicate the presence of clumps 
at a distance of 0.1 pc. The flat central region in the profile cor- 
responds to the Jeans length below which the effect of gravity is 
negligible. The overall density profiles agree on the scale above 0. 1 
pc although there are differences in the interior, which are likely 
related to the morphological differences discussed below. 

The upper right panel of figure [T] depicts the temperature ra- 
dial profile. The temperature initially rises at larger radii due to the 
virialization shocks and reaches a maximum around 10 4 K which 
corresponds to the virial temperature of the halo. The cooling due 
to Lyman alpha radiation becomes important at these temperatures 
and brings the temperature of the halo down to about 7000 K. In the 
presence of strong UV Lyman Werner flux J 2 | = 10 3 the formation 
of H 2 remains inhibited and the halo collapses almost isothermally. 
It is found that the thermal properties of the halo are independent 
of the Jeans resolution. The total energy radial profile is shown in 
the bottom left panel of figureQ] The value of total energy is a few 
times 10 12 erg/g which is comparable to our previous study. There 
is an increase in the total energy at outer radii due to the infall of 
the gas. It becomes almost constant at intermediate radii. The vari- 
ations in the core may arise from the differences in the morphology 



4 M. A. Latif et al. 




Figure 2. The figure illustrates the state of simulations at the collapse redshift. Density projections are shown for the central 4000 AU region of the halo 
for various Jeans resolutions. The top panels show the cases for 16 (left) and 32 (right) cells per Jeans length. The bottom panels depict the cases for Jeans 
resolutions of 64 and 128 cells from left to right respectively. 



of the halo. On scales above 0. 1 pc, the total energy is converged 
for different resolutions. 

The averaged radial profile for the magnetic field strength is 
presented in the bottom right panel of the figure [TJ The simula- 
tions were seeded with the same initial magnetic field strength of 
1CT 14 G at the initial redshift. The values of the magnetic field fol- 
low a power law behavior for different resolutions at larger radii. 
The magnetic field amplification seems to proceed via gravitational 
compression at larger radii while turbulent amplification kicks in in 
the interior. The maximum value of the magnetic field strength is 
about 10~ 4 G. The value of B raae is about two orders of magnitude 
higher for the J128 case compared to the Jie- This increase in mag- 
netic field strength is a consequence of the turbulent amplification 
of the magnetic field as turbulent eddies are better resolved in the 
highest resolution case. We will focus on the details of amplifica- 
tion in the next subsection. 



3.1.2 Morphology of the halo 

The state of simulations at the collapse redshift for different Jeans 
resolutions is illustrated in the density projections of figure [2] Re- 
markable dissimilarities in the morphology of the halo are found 
for different resolutions. For the J 16 case, the halo is collapsed into 
a single filament. There are no turbulent structures present in this 
case which indicates that turbulence is clearly under-resolved. The 



more complex clumps are observed for a Jeans resolution of > 
32 cells. They are very well separated from each other and may 
form binary or even multiple systems. Again, turbulence seems to 
be marginally resolved in this case. For the Jeans resolution of 64 
cells, the medium inside the halo becomes turbulent. One can see 
that again there are multiple clumps and turbulence is better re- 
solved compared to the previous cases. They are not gravitationally 
bound and have sub-solar masses. 



For the highest resolution case, turbulence seems to be very 
well resolved. There are no coherent structures present in the halo. 
Our results suggest that one needs Jeans resolution of at least > 64 
cells to resolve gravity driven turbulence in atomic cooling haloes. 
In the highest resolution study case, the halo has become highly 
turbulent. This can have important implications for the ultimate 
fate of the halo. Turbulence can efficiently transport angular mo- 
mentum and may l ead to the formation of a see d black hole or a 
supermassive star jBegelman & Shlosman|[2009h . Due to exorbi- 
tant computational costs, we stopped simulations soon after they 
reached the maximum refinement as further evolution would pro- 
vide unphysical results. Nonetheless, our results show structure of 
the halo strongly varies with resolution and the presence of turbu- 
lence may lead to different fragmentation properties. 
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Figure 3. Density-weighted magnetic field strength projections for various Jeans resolutions are shown in this figure. They are depicted at the central field of 
view of 4000 AU. Top panels show the cases for 16 (left) and 32 (right) cells per Jeans length. Bottom panels depict the cases for Jeans resolutions of 64 and 
128 cells from left to right respectively. 
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Figure 4. In the left panel of the figure, magnetic energy scaled by p 4 ' 3 is plotted against the gas density for various Jeans resolutions as indicated in the 
legend. Radial profiles of the vorticity normalized by the free-fall time are shown in the right panel. 
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Figure 5. This figure depicts the amplification of magnetic field strength by shear (left) and compression (right). The solid lines represent the contribution 
of positive components while dotted lines stand for the negative components of both shear and compression. Blue, green, purple and red colors represent the 
Jeans resolutions of 16, 32, 64 and 128 cells respectively. 
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Figure 6. In the left panel of the figure turbulent Mach number is plotted against gas density. The radial profile of turbulent viscosity is shown in the right 
panel. Different colors represent various resolutions as indicated in legends. 



3.1.3 Magnetic field Amplification and turbulence 

The simulations were started with the same initial seed magnetic 
field of 10~ 14 G (proper) at redshift 99. The magnetic field is ampli- 
fied during the collapse, reaching a peak value of 100 G at densities 
of 10~ 12 g/cm 3 at z = 1 1.5. In total, the magnetic field is amplified 
about 10 orders of magnitude. This is evident in figure[3]where the 
total strength of the magnetic field is depicted for different Jeans 
resolutions. The strength of magnetic field is correlated with the 
density structure present inside the halo. The amplification of mag- 
netic fields for Ji6 and J32 cases is mainly due to the gravitational 
compression as amount of turbulence is poorly resolved. 

In fact, we clearly recognize turbulent structures for Jeans res- 
olutions of > 64 cells. The turbulent energy is converted into mag- 
netic energy via the small scale dynamo and exponentially ampli- 
fies the magnetic field in addition to the maximum possible con- 
tribution from flux freezing. The large region with higher values 
of vorticity in the J128 case enhances the magnetic energy which 



results in enhanced magnetic field strengths as shown in the bot- 
tom right panel of figure[3] The magnetic field strength is higher in 
the center as turbulence is produced locally due to the gravitational 
collapse. Our results clearly demonstrate the necessity of a Jeans 
length resolution of > 64 cells to excite the dynamo action. The 
critical resolution is thus higher than in previous studies, wh ich is 
likely due to the higher Mach numbers jFederrath et alj|201 ll) . 

We have further quantified the contribution of the dynamo 
generated magnetic field as manifested in the left panel of figure 
[4] The magnetic energy would scale as E B = p 4/3 according to 
the expectation of spherical collapse. We infer the dynamo ampli- 
fication of the magnetic field by dividing the magnetic energy by 
the maximum possible contribution from the frozen magnetic field 
lines. It can be noticed from the plot that there is no contribution of 
dynamo action for a J, 6 run and magnetic energy seems to be dissi- 
pated. For higher resolution runs, initially there is a slight decrease 
in magnetic energy, which can be due to the strong virialization 



Magnetic field amplification in massive primordial haloes 7 



shocks that destroy the coherence in turbulence and magnetic fields 
necess a ry to derive the d ynamo (as interpreted bv lFederrath et al.l 
feOllh : ISur et al.l \lQ\± \ Later on the contribution of the small 
scale dynamo becomes important and efficiently amplifies the mag- 
netic field. For the best resolution case (Ji2s)> dynamo amplification 
is about two orders of magnitude above the flux freezing. 

In the right panel of figure [6] the Mach number is plotted 
against the gas density for different Jeans resolutions. The Mach 
number peaks around three. For the highest resolution case, the 
maximum Mach number in the core of the halo is two. Over- 
all, turbulence inside the halo is supersonic. The small scale dy- 
namo is known to be less efficient for c ompressible turbulence 
jFederrath et al ] |201lUSchoberetalJl2012h . as observed in the ha- 
los modeled here. The amplification of the magnetic field is there- 
fore less ef ficient compared to the minihaloes explored in previ- 
ous studies dSur et al.l2O10l ; iFederrath et alj201 ll ; iTurk et alj20H 
ISchober et alj|2012f) . However, the latter is likely compensated in 
astrophysical systems by the much hi gher Reynolds numbe rs, as 
the growth rate scales as Re 1 ' 3 - Re 1/2 dSchober et al 1l2012al) . 

Additional evidence for the turbulent amplification of mag- 
netic fields is shown in the right panel of figure|4] The radial profile 
of vorticity normalized by the local free-fall time demonstrates that 
the amount of vorticity increases with resolution particularly in the 
center of the halo. The relatively high vorticity in the inner region 
indicates turbulence driven by the collapse. The peak in vorticity at 
larger radii is because of virialization shocks. 

Gravitational compression squeezes magnetic field lines un- 
der the condition of flux freezing for the ideal MHD conditions. 
The rate of change of magnetic pressure can be computed from the 
induction equation (see Schmidt et al 2012) 

IO=i( B ' B ^- 2/3B2d )' a) 

where §- is 4- + v.V and S* = Sh _ 1/35;; is the trace-free rate 

Dt at ij J J 

of strain tensor. The first term on the right hand side of equation [T] 
comes from the action of turbulent shear on magnetic fields while 
the second term results from the gas compression both due to grav- 
ity and shocks. In each case the positive contribution exceeds the 
negative contribution, implying a net amplification effect. We have 
computed the contributions of both shear and compression for dif- 
ferent resolutions and they are plotted in figure [5] We see a clear 
trend of enhanced amplification rate with increasing resolution. The 
amplification by local shear and compression is quite comparable. 
This has already been seen in MHD turbulence simulations and 
seems to be a general property of self-gravitating turbulence. More- 
over, amplification correlates very well with density as expected in 
gravity driven turbulence of collapsing halos (Schmidt et al 2012). 

We have also computed the turbulent viscosity ( for definition 
see equation [5} for different Jeans resolutions, which we estimate 
as the product of Jeans length and the turbulent velocity. We expect 
that to be a reasonable approximation, as turbulence is produced 
as a result of gravitational collapse. The turbulent velocity is com- 
puted as : 

V turb = ^'lu ~ V l-, d - V ml (2) 

where v rms is the root mean square velocity, v rat j is the radial ve- 
locity and v rot is the rotational velocity. The radial profile of the 
turbulent viscosity is depicted in figure [6] showing a power law 
behavior regardless of the resolution and decreases towards the 
smaller radii. Although one would expect the turbulent viscosity to 
decrease with increasing resolution, the behavior is more compli- 



cated. The amount of viscosity is enhanced in the center for higher 
resolution runs because of better resolved turbulent structures. The 
maximum value of the turbulent viscosity is few times 10 9 cm 2 /s 
for the best resolution case. 

3.2 Saturation of magnetic field 

Only a modest range of Reynolds numbers can be achieved in nu- 
merical simulations due to the exorbitant computational costs. Sim- 
ulations starting from weak fields thus will not achi eve saturation, 
as they will still underestimate the growth rate (see ISchober et al.l 
d2012bh ). Given the constraints in achieving saturation of magnetic 
fields in numerical experiments, we conducted a number of simula- 
tions with different initial fields as listed in table [T] The properties 
of the halo for different initial seed fields are shown in the figure[7J 
Similar to the previous cases, the density follows a power law be- 
havior (i.e., R~ 2 ) and the halo collapses isothermally as seen in the 
temperature radial profile. The total energy profile is analogous to 
the earlier cases as depicted in the figure|7j The variations in the in- 
terior of the halo arise due to the dissimilarities in the morphology 
for different field strengths. The strength of the magnetic field for 
various seed fields is depicted in the bottom right panel. Magnetic 
field is amplified up to ~ 0. 1 G for the highest seed field. 

The morphology of the halo is strongly influenced by the 
strength of the magnetic field as demonstrated in figure [8] In gen- 
eral, it is found that the enhanced strength of the magnetic field 
leads to the formation of compact filamentary structures and seems 
to depreciate the role of turbulence. Perhaps, it suggests that mag- 
netic field is close to the saturation. 

To assess the saturation of the dynamo generated magnetic 
field in a more elaborate way, we have plotted the ratio of magnetic 
to kinetic energy for these runs which is shown in the left panel 
of figure|9] This ratio increases with density for all seed fields. We 
infer the saturation of magnetic field from the change in slope of 
Eb/Ek, as the latter will change once dynamo amplification stops. 
The magnetic energy associated with small scale dynamo attains 
the peak value about 5 x 10~ 3 the total kinetic energy and then starts 
to decrease. At that point, we reach approximate equipartition with 
the turbulent energy. This behavior is obtained starting from a seed 
field of 5x 10" 11 G. To elucidate the saturation of dynamo amplified 
magnetic field, we have plotted the ratio of magnetic to turbulent 
energy and is depicted in the right panel of figure [9] The turbulent 
energy is given by 

1 , 

E.urb = -pV-, rb (3) 

here p is the gas mass density and v lllr b is the turbulent velocity. 
The plot clearly demonstrates that magnetic energy comes close to 
equipartition with turbulent energy and dynamo amplification gets 
saturated. Nevertheless, the amplification by gravitational com- 
pression will continue until back reactions occur and then system 
evolve towards a complete saturati on. In general, our results are in 
agree ment with previous studies dFederrath et alj|201 ll ; ISur et al.l 
|2012|) . They also find similar values for the saturation of magnetic 
energy. The saturation value depends on the turbulent energy injec- 
tion scale, type of driving and the Mach number of turbulence. 

4 DISCUSSION AND CONCLUSIONS 

In all, we have conducted 1 1 cosmological magnetohydrodynam- 
ical simulations using the adaptive mesh refinement code Enzo to 
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Figure 7. This figure shows the radially binned spherically averaged radial profiles for various initial seed fields as shown in the legend. The upper left panel 
of the figure shows the density radial profile. The temperature radial profile is depicted the upper right panel. The bottom left panel of the figure shows the 
total energy radial profile. The radial profile for the total magnetic field strength is shown in the bottom right panel. 



Table 1. The list of simulations performed to probe the saturation of mag- 
netic fields 



Model Magnetic field strength 
[Gauss] 



Bl 


1 X IO-' 4 


B2 


6x 1CT 12 


B? 


1.8 X IO"" 


B4 


5.4 X IO"" 


B5 


6.7 X 10-" 


B6 


2x 10 -10 



study the magnetization of massive primordial haloes. We started 
our simulations with nested grid initial conditions and inserted 27 
additional refinement levels during the course of simulations. It 
gives us an effective resolution down to sub AU scales. To explore 
the turbulent amplification of magnetic fields in detail, we man- 
dated a Jeans refinement of 16, 32, 64 and 128 cells. These simu- 
lations were started at redshift 99 with initial seed field of 10" 14 G. 
We also probed the saturation of magnetic field for a fixed Jeans 
resolution of 64 cells. Here we briefly summarize our main find- 
ings: 



• No turbulent amplification on large scales due to strong viri- 
alization shocks. 

• The turbulent dynamo becomes active in the interior of the 
halo and amplifies the magnetic field. 

• The amplification efficiency is reduced compared to the mini- 
haloes due to the higher Mach numbers, but likely compensated by 
the Reynolds numbers in real system. 

• The saturation of magnetic field occurs for the initial seed field 
of ~ 5 x 10-" G. 

• The morphology of the halo is influenced when the magnetic 
field is close to the saturation state. 

In all simulations, the gas falls into the dark matter potential 
and gets heated up to its virial temperature. At temperatures > 10 4 
K, the cooling due to Lyman alpha comes into play and brings the 
temperature of the gas down to 7000 K. The formation of H2 re- 
mains suppressed due to the presence of strong photo-dissociation 
UV radiation flux (i.e., J 2 i = 10 3 ). Consequently, almost isothermal 
collapse occurs. We reach a maximum peak density of 10-' 2 g/cm 3 
in our simulation. The density radial profile follows an isothermal 
behavior R 2 according to the expectation of isoth ermal collapse. 
These results are in agreement with previous studies dOh & Haimanl 
l2002MRegan & Haehneltl2009cMWise et alj20O8tlLatif et al.l2012l) 
and suggest that thermal properties of the haloes are independent of 
resolution effects and magnetic seed fields. 
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Figure 8. Figure shows the density projections at the collapse redshift for various initial seed fields. The central 10000 AU region of the halo is shown for a 
fixed resolution of 64 cells per Jeans length. The top panels represent Bl (left) and B2 (right), middle panels B3 (left) and B4 (right), and bottom panels B5 
(left) and B6 (right). 
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We found that the average properties of the halo are almost 
converged for different resolutions whereas strong dissimilarities in 
the morphology of the halo are observed for different Jeans resolu- 
tion. The turbulence is very well resolved for the best resolution run 
(i.e., J 12 g) and will have potential implication for the birth of black 
holes via direct collapse scenario. Our study suggest higher values 
of Jeans refinement for atomic cooling haloes (i.e., > J64) because 
of higher Mach numbers and larger spatial scales. Previous stud- 
ies found lower Jeans resolution > 32 cells for minihaloes where 
turbu l ence is typically subsonic dLatifetalj|2012l ;l Federrath et al .1 
1201 ll : iTurk et alj|2012h . Turbulence causes the medium inside the 
halo to fragment into clumps and this process is very sensitive to 
resolution. We stopped simulations because of computational con- 
straints. In future studies techniques like Jeans heating or sink par- 
ticles would be employed to advance simulation for longer dynam- 
ical times scales. 

We examined the amplification of magnetic fields and found 
that it is amplified about ten orders of magnitude. The maximum 
field strength (i.e., J 12 8) is about 100 /j Gauss. It is initially ampli- 
fied by the gravitational compression under the condition of flux 
freezing for low resolution runs (i.e, J 16 and J 32 ). The role of tur- 
bulent amplification at radii larger than 0. 1 pc is depreciated due to 
strong virialization shocks. Turbulence is mainly generated during 
the gravitational collapse. It drives the small scale dynamo which 
consequently amplifies the magnetic field strength by randomly 
stretching, twisting and folding magnetic lines. In the kinematic 
regime one expects that the growth of magnetic energy is exponen- 
tial in time while it becomes linear once backreactions of the mag- 
netic field start playing a role. For the highest resolution case (i.e., 
Ji2g), the magnetic field is amplified by at least an order of magni- 
tude above the maximum possible amplification by flux freezing. 
Indeed, the small scale dynamo gets excited by the gravitational 
collapse of cosmological initial conditions provided that resolution 
is enough to capture turbulent vortical motions. In general, ampli- 
fication of the magnetic field is less efficient compared to the mini- 
haloes because of the higher Mach numbers. Astrophysical systems 
however have higher Reynolds numbers than those realized in nu- 
merical simulations. As the growth rate scales as Re 1 ' 3 - Re 1/2 , the 



amplifi cation found here is thus only a lower limit l lSchober et al .1 
l2012ah . 

We explored the expected saturation values of magnetic fields 
in MHD simulations by selecting initial seed fields of various 
strengths. The ratio of magnetic to kinetic energy reaches a peak 
value of 5 x 10~ 3 times equipartition value and then starts to de- 
crease when the magnetic energy becomes comparable to the tur- 
bulent energy. Further amplification via compression will however 
occur during gravitational collapse. We note here that the enhanced 
field strength also influences the morphology of the halo, as the 
structure becomes compact for stronger seed fields. 

The simulations above confirm that the small-scale dynamo is 
operational in massive primordial halos, which are generally con- 
sidered as candidates for black hole formation via direct collapse. 
During the kinematic phase, i.e. before backreactions occur, a cen- 
tral property of the dynamo is an exponential growth, with a rate 
that scales as Re 112 for Burgers tu rbulence l lSubramanianl 1 19991 : 
iBrandenburg & Subramanianl2005l) . and Re 1/3 in the extreme case 
of Kolmogorov turbulence. As we are far from resolving the vis- 
cous scales, one can thus expect growth rates considerably larger 
than those obtained in numerical simulations. Consequently, the 
smallest scales saturate first, and the dynamo growth proceeds on 
larger scales, for which the magnetic field strength is not sufficient 
to suppress further amplification. In this non-linear regime, satura- 
tion is expected within a few to a few doz en eddy- turnover times 
dSchekochihin et al.ll2002klBeresnvakll2012h . 

Thus, we expect the presence of magnetic fields in approx- 
imate equipartition with turbulent energy in the centers of such 
halos. The latter may aid the formation and subsequent growth 
of massive black holes both via their magnetic pressure and the 
magnetic tension, providing an additional means for angular mo- 
mentum transport in the accretion disk. While the magnetic pres- 
sure is best described via a magnetic Jeans mass, given as Mj B oc 
B 3 /p 2 , the effect of magnetic tension is likely more prominent af- 
ter the formation of accre tion disks. Adoptin g the framework of 
Shakura & Sunvaevl i 19731) and lPringld dl98lh . we describe the ac- 
cretion rate in such a disk via 
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M(r) = InvLJr) 



rflnn(r) 



cl\nr 



(4) 



with S e (r) the gas surface density of the disk, fi(r) the angular ve- 
locity, and v an effective viscosity, given as 



v = av, urb (r)h(r). 



(5) 



Here, v turb denotes the turbulent velocity, and h(r) the scale-height 
of the disk. The parameter a denotes the efficiency of the accretion 
process, and has contributions from turbule nce, tangled magnetic 
fields and gravitational stresses. It is given as dBalbus & Papaloizoul 
1 19991; lHawlevll2000l : iFromang et al1l2004l) 



y,Q + 7ff w V,Q + r%>,0 

P(r,t) 



(6) 



where P denotes the thermal pressure of the gas and T r ^, represents 
the different contributions of stresses. While initially, the contri- 
bution of turbulence and magnetic fields to the angular momen- 
tum transport is likely comparable, we expect the magnetic field 
to be amplified further in the accretion disk, as a result of the 
a — SI dynamo dBrandenburg & Subramanianll2005h . As discussed 
by IPudritzl d 198 lh . the build-up of a larger-scale magnetic field 
component is likely accompanied by an increase in the fluctuations 
of B. The presenc e of a strong coherent field will further give rise to 
jets a nd outflows dCamenzin d 1995; Mc Kinney 2005; Pudritz et all 
|2012|) . These provide a highly efficient mechanism of transporting 
angular momentum. Considering that the rotation timescale of such 
accretion disks is considerably smaller than the Hubble time, these 
disks are likely to reach strong coherent fields well before redshift 
6. The latter explains the strong rotation measures observed for 
high-redshift quasars, and potentially enhances the accretion rate 
to form more massive black holes. In the central accretion disk, 
coherent magnetic field structures are likely produced well before 
redshift 6, implying the presence of strong rotation measures and 
enhanced accretion rates. 
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